Multi-omic Analysis of Uterine Leiomyomas in Self-Described Black and White Women: Molecular Insights into Health Disparities.

BACKGROUND: Black women are at an increased risk to develop uterine leiomyomas (ULMs) and to experience worse disease prognosis compared to White women. Epidemiological and molecular factors have been identified as underlying these disparities, but there remains a paucity of deep, multi-omic analysis investigating molecular differences in ULMs from Black and White patients. OBJECTIVE: To identify molecular alterations within ULM tissues correlating with patient race by multi-omic analyses of ULMs collected from cohorts of Black and White women. STUDY DESIGN: We performed multi-omic analysis of ULMs from Black (42) and White (47) women undergoing hysterectomy for symptomatic uterine leiomyomata. Our analysis also included application of orthogonal methods to evaluate fibroid biomechanical properties, such as second harmonic generation microscopy, uniaxial compression testing, as well as sheer wave ultrasonography analyses. RESULTS: We found a greater proportion of mediator complex subunit 12 (MED12) mutant ULMs from Black women (>35% increase, Mann Whitney U p = 7E-4). MED12 mutant tumors exhibited elevated abundance of extracellular matrix proteins, including several collagen isoforms, involved in regulation of the core matrisome. Histological analysis of tissue fibrosis using trichrome staining and secondary harmonic generation microscopy confirmed that MED12 mutant tumors are more fibrotic than MED12 wildtype tumors. Using Sheer Wave ultrasonography in a prospectively collected cohort, Black patients had fibroids that were firmer when compared to White patients, even when similar in size. These analyses also uncovered expression quantitative trait loci (eQTL) associated with altered allele frequencies in African and European populations that correlated with differential abundance of several proteins in ULM that are independent of MED12 mutational status, including multiple eQTL mapping to tetratricopeptide repeat protein 38. CONCLUSIONS: Our study shows that Black women have a higher prevalence of ULMs harboring mutations in MED12 and that this mutational status correlates with increased tissue fibrosis in comparison with wildtype ULMs. Our study provides insights into molecular alterations underlying racial disparities in ULMs and improves our understanding of the molecular etiology underlying ULM development within these populations.

Immune response profiles from humans experimentally exposed to Phlebotomus duboscqi bites.

Introduction: Cutaneous leishmaniasis is a neglected vector-borne parasitic disease prevalent in 92 countries with approximately one million new infections annually. Interactions between vector saliva and the human host alter the response to infection and outcome of disease. Methods: To characterize the human immunological responses developed against saliva of Phlebotomus duboscqi, a Leishmania major (L. major) vector, we repeatedly exposed the arms of 14 healthy U.S volunteers to uninfected P. duboscqi bites. Blood was collected a week after each exposure and used to assess total IgG antibodies against the proteins of P. duboscqi salivary gland homogenate (SGH) and the levels of IFN-gamma and IL-10 from peripheral blood mononuclear cells (PBMCs) stimulated with SGH or recombinant sand fly proteins. We analyzed skin punch biopsies of the human volunteer arms from the insect bite site and control skin site after multiple P. duboscqi exposures (four volunteers) using immunohistochemical staining. Results: A variety of immediate insect bite skin reactions were observed. Late skin reactions to insect bites were characterized by macular hyperpigmentation and/or erythematous papules. Hematoxylin and eosin staining showed moderate mononuclear skin infiltrate with eosinophils in those challenged recently (within 2 months), eosinophils were not seen in biopsies with recall challenge (6 month post bites). An increase in plasma antigen-specific IgG responses to SGH was observed over time. Western Blot results showed strong plasma reactivity to five P. duboscqi salivary proteins. Importantly, volunteers developed a cellular immunity characterized by the secretion of IFN-gamma upon PBMC stimulation with P. duboscqi SGH and recombinant antigens. Discussion: Our results demonstrate that humans mounted a local and systemic immune response against P. duboscqi salivary proteins. Specifically, PduM02/SP15-like and PduM73/adenosine deaminase recombinant salivary proteins triggered a Th1 type immune response that might be considered in future development of a potential Leishmania vaccine.

Identification of germline population variants misclassified as cancer-associated somatic variants.

Introduction: Databases used for clinical interpretation in oncology rely on genetic data derived primarily from patients of European ancestry, leading to biases in cancer genetics research and clinical practice. One practical issue that arises in this context is the potential misclassification of multi-ancestral population variants as tumor-associated because they are not represented in reference genomes against which tumor sequencing data is aligned. Methods: To systematically find misclassified variants, we compared somatic variants in census genes from the Catalogue of Somatic Mutations in Cancer (COSMIC) V99 with multi-ancestral population variants from the Genome Aggregation Databases' Linkage Disequilibrium (GnomAD). By comparing genomic coordinates, reference, and alternate alleles, we could identify misclassified variants in genes associated with cancer. Results: We found 192 of 208 genes in COSMIC's cancer-associated census genes (92.31%) to be associated with variant misclassifications. Among the 1,906,732 variants in COSMIC, 6,957 variants (0.36%) aligned with normal population variants in GnomAD, concerning for misclassification. The African / African American ancestral population included the greatest number of misclassified variants and also had the greatest number of unique misclassified variants. Conclusion: The direct, systematic comparison of variants from COSMIC for co-occurrence in GnomAD supports a more accurate interpretation of tumor sequencing data and reduces bias related to genomic ancestry.

ProteoMixture: A cell type deconvolution tool for bulk tissue proteomic data.

Numerous multi-omic investigations of cancer tissue have documented varying and poor pairwise transcript:protein quantitative correlations, and most deconvolution tools aiming to predict cell type proportions (cell admixture) have been developed and credentialed using transcript-level data alone. To estimate cell admixture using protein abundance data, we analyzed proteome and transcriptome data generated from contrived admixtures of tumor, stroma, and immune cell models or those selectively harvested from the tissue microenvironment by laser microdissection from high grade serous ovarian cancer (HGSOC) tumors. Co-quantified transcripts and proteins performed similarly to estimate stroma and immune cell admixture (r ≥ 0.63) in two commonly used deconvolution algorithms, ESTIMATE or ConsensusTME. We further developed and optimized protein-based signatures estimating cell admixture proportions and benchmarked these using bulk tumor proteomic data from over 150 patients with HGSOC. The optimized protein signatures supporting cell type proportion estimates from bulk tissue proteomic data are available at https://lmdomics.org/ProteoMixture/.

Proteogenomic analysis of enriched HGSOC tumor epithelium identifies prognostic signatures and therapeutic vulnerabilities.

We performed a deep proteogenomic analysis of bulk tumor and laser microdissection enriched tumor cell populations from high-grade serous ovarian cancer (HGSOC) tissue specimens spanning a broad spectrum of purity. We identified patients with longer progression-free survival had increased immune-related signatures and validated proteins correlating with tumor-infiltrating lymphocytes in 65 tumors from an independent cohort of HGSOC patients, as well as with overall survival in an additional 126 HGSOC patient cohort. We identified that homologous recombination deficient (HRD) tumors are enriched in pathways associated with metabolism and oxidative phosphorylation that we validated in independent patient cohorts. We further identified that polycomb complex protein BMI-1 is elevated in HR proficient (HRP) tumors, that elevated BMI-1 correlates with poor overall survival in HRP but not HRD HGSOC patients, and that HRP HGSOC cells are uniquely sensitive to BMI-1 inhibition.

Metronomic dosing of ovarian cancer cells with the ATR inhibitor AZD6738 leads to loss of CDC25A expression and resistance to ATRi treatment.

OBJECTIVE: ATR kinase inhibitors promote cell killing by inducing replication stress and through potentiation of genotoxic agents in gynecologic cancer cells. To explore mechanisms of acquired resistance to ATRi in ovarian cancer, we characterized ATRi-resistant ovarian cancer cells generated by metronomic dosing with the clinical ATR inhibitor AZD6738. METHODS: ATRi-resistant ovarian cancer cells (OVCAR3 and OV90) were generated by dosing with AZD6738 and assessed for sensitivity to Chk1i (LY2603618), PARPi (Olaparib) and combination with cisplatin or a CDK4/6 inhibitor (Palbociclib). Models were characterized by diverse methods including silencing CDC25A in OV90 cells and assessing impact on ATRi response. Serum proteomic analysis of ATRi-resistant OV90 xenografts was performed to identify circulating biomarker candidates of ATRi-resistance. RESULTS: AZD6738-resistant cell lines are refractory to LY2603618, but not to Olaparib or combinations with cisplatin. Cell cycle analyses showed ATRi-resistant cells exhibit G1/S arrest following AZD6738 treatment. Accordingly, combination with Palbociclib confers resistance to AZD6738. AZD6738-resistant cells exhibit altered abundances of G1/S phase regulatory proteins, including loss of CDC25A in AZD6738-resistant OV90 cells. Silencing of CDC25A in OV90 cells confers resistance to AZD6738. Serum proteomics from AZD6738-resistant OV90 xenografts identified Vitamin D-Binding Protein (GC), Apolipoprotein E (APOE) and A1 (APOA1) as significantly elevated in AZD6738-resistant backgrounds. CONCLUSIONS: We show that metronomic dosing of ovarian cancer cells with AZD6738 results in resistance to ATR/ Chk1 inhibitors, that loss of CDC25A expression represents a mechanism of resistance to ATRi treatment in ovarian cancer cells and identify several circulating biomarker candidates of CDC25A low, AZD6738-resistant ovarian cancer cells.

NSC243928 Treatment Induces Anti-Tumor Immune Response in Mouse Mammary Tumor Models.

NSC243928 induces cell death in triple-negative breast cancer cells in a LY6K-dependent manner. NSC243928 has been reported as an anti-cancer agent in the NCI small molecule library. The molecular mechanism of NSC243928 as an anti-cancer agent in the treatment of tumor growth in the syngeneic mouse model has not been established. With the success of immunotherapies, novel anti-cancer drugs that may elicit an anti-tumor immune response are of high interest in the development of novel drugs to treat solid cancer. Thus, we focused on studying whether NSC243928 may elicit an anti-tumor immune response in the in vivo mammary tumor models of 4T1 and E0771. We observed that NSC243928 induced immunogenic cell death in 4T1 and E0771 cells. Furthermore, NSC243928 mounted an anti-tumor immune response by increasing immune cells such as patrolling monocytes, NKT cells, B1 cells, and decreasing PMN MDSCs in vivo. Further studies are required to understand the exact mechanism of NSC243928 action in inducing an anti-tumor immune response in vivo, which can be used to determine a molecular signature associated with NSC243928 efficacy. NSC243928 may be a good target for future immuno-oncology drug development for breast cancer.

Lymphocyte antigen 6K signaling to aurora kinase promotes advancement of the cell cycle and the growth of cancer cells, which is inhibited by LY6K-NSC243928 interaction.

Lymphocyte antigen 6K (LY6K) is a small GPI-linked protein that is normally expressed in testes. Increased expression of LY6K is significantly associated with poor survival outcomes in many solid cancers, including cancers of the breast, ovary, gastrointestinal tract, head and neck, brain, bladder, and lung. LY6K is required for ERK-AKT and TGF-ß pathways in cancer cells and is required for in vivo tumor growth. In this report, we describe a novel role for LY6K in mitosis and cytokinesis through aurora B kinase and its substrate histone H3 signaling axis. Further, we describe the structural basis of the molecular interaction of small molecule NSC243928 with LY6K protein and the disruption of LY6K-aurora B signaling in cell cycle progression due to LY6K-NSC243928 interaction. Overall, disruption of LY6K function via NSC243928 led to failed cytokinesis, multinucleated cells, DNA damage, senescence, and apoptosis of cancer cells. LY6K is not required for vital organ function, thus inhibition of LY6K signaling is an ideal therapeutic approach for hard-to-treat cancers that lack targeted therapy such as triple-negative breast cancer.

Integrated multi-omic analysis of low-grade ovarian serous carcinoma collected from short and long-term survivors.

BACKGROUND: Low-grade serous ovarian cancer (LGSOC) is a rare disease that occurs more frequently in younger women than those with high-grade disease. The current treatment is suboptimal and a better understanding of the molecular pathogenesis of this disease is required. In this study, we compared the proteogenomic analyses of LGSOCs from short- and long-term survivors (defined as < 40 and > 60 months, respectively). Our goal was to identify novel mutations, proteins, and mRNA transcripts that are dysregulated in LGSOC, particularly in short-term survivors. METHODS: Initially, targeted sequencing of 409 cancer-related genes was performed on 22 LGSOC and 6 serous borderline ovarian tumor samples. Subsequently, whole-genome sequencing analysis was performed on 14 LGSOC samples (7 long-term survivors and 7 short-term survivors) with matched normal tissue samples. RNA sequencing (RNA-seq), quantitative proteomics, and phosphoproteomic analyses were also performed. RESULTS: We identified single-nucleotide variants (SNVs) (range: 5688-14,833 per sample), insertion and deletion variants (indels) (range: 880-1065), and regions with copy number variants (CNVs) (range: 62-335) among the 14 LGSOC samples. Among all SNVs and indels, 2637 mutation sites were found in the exonic regions. The allele frequencies of the detected variants were low (median12%). The identified recurrent nonsynonymous missense mutations included KRAS, NRAS, EIF1AX, UBR5, and DNM3 mutations. Mutations in DNM3 and UBR5 have not previously been reported in LGSOC. For the two samples, somatic DNM3 nonsynonymous missense mutations in the exonic region were validated using Sanger sequencing. The third sample contained two missense mutations in the intronic region of DNM3, leading to a frameshift mutation detected in RNA transcripts in the RNA-seq data. Among the 14 LGSOC samples, 7754 proteins and 9733 phosphosites were detected by global proteomic analysis. Some of these proteins and signaling pathways, such as BST1, TBXAS1, MPEG1, HBA1, and phosphorylated ASAP1, are potential therapeutic targets. CONCLUSIONS: This is the first study to use whole-genome sequencing to detect somatic mutations in LGSOCs with matched normal tissues. We detected and validated novel mutations in DNM3, which were present in 3 of the 14 samples analyzed. Additionally, we identified novel indels, regions with CNVs, dysregulated mRNA, dysregulated proteins, and phosphosites that are more prevalent in short-term survivors. This integrated proteogenomic analysis can guide research into the pathogenesis and treatment of LGSOC.

Proteogenomic analysis of lung adenocarcinoma reveals tumor heterogeneity, survival determinants, and therapeutically relevant pathways.

We present a deep proteogenomic profiling study of 87 lung adenocarcinoma (LUAD) tumors from the United States, integrating whole-genome sequencing, transcriptome sequencing, proteomics and phosphoproteomics by mass spectrometry, and reverse-phase protein arrays. We identify three subtypes from somatic genome signature analysis, including a transition-high subtype enriched with never smokers, a transversion-high subtype enriched with current smokers, and a structurally altered subtype enriched with former smokers, TP53 alterations, and genome-wide structural alterations. We show that within-tumor correlations of RNA and protein expression associate with tumor purity and immune cell profiles. We detect and independently validate expression signatures of RNA and protein that predict patient survival. Additionally, among co-measured genes, we found that protein expression is more often associated with patient survival than RNA. Finally, integrative analysis characterizes three expression subtypes with divergent mutations, proteomic regulatory networks, and therapeutic vulnerabilities. This proteogenomic characterization provides a foundation for molecularly informed medicine in LUAD.

Whole-Genome Sequencing Identifies PPARGC1A as a Putative Modifier of Cancer Risk in BRCA1/2 Mutation Carriers.

While BRCA1 and BRCA2 mutations are known to confer the largest risk of breast cancer and ovarian cancer, the incomplete penetrance of the mutations and the substantial variability in age at cancer onset among carriers suggest additional factors modifying the risk of cancer in BRCA1/2 mutation carriers. To identify genetic modifiers of BRCA1/2, we carried out a whole-genome sequencing study of 66 ovarian cancer patients that were enriched with BRCA carriers, followed by validation using data from the Pan-Cancer Analysis of Whole Genomes Consortium. We found PPARGC1A, a master regulator of mitochondrial biogenesis and function, to be highly mutated in BRCA carriers, and patients with both PPARGC1A and BRCA1/2 mutations were diagnosed with breast or ovarian cancer at significantly younger ages, while the mutation status of each gene alone did not significantly associate with age of onset. Our study suggests PPARGC1A as a possible BRCA modifier gene. Upon further validation, this finding can help improve cancer risk prediction and provide personalized preventive care for BRCA carriers.

Germline mutation landscape of DNA damage repair genes in African Americans with prostate cancer highlights potentially targetable RAD genes.

In prostate cancer, emerging data highlight the role of DNA damage repair genes (DDRGs) in aggressive forms of the disease. However, DDRG mutations in African American men are not yet fully defined. Here, we profile germline mutations in all known DDRGs (N = 276) using whole genome sequences from blood DNA of a matched cohort of patients with primary prostate cancer comprising of 300 African American and 300 European Ancestry prostate cancer patients, to determine whether the mutation status can enhance patient stratification for specific targeted therapies. Here, we show that only 13 of the 46 DDRGs identified with pathogenic/likely pathogenic mutations are present in both African American and European ancestry patients. Importantly, RAD family genes (RAD51, RAD54L, RAD54B), which are potentially targetable, as well as PMS2 and BRCA1, are among the most frequently mutated DDRGs in African American, but not in European Ancestry patients.

Peptide ancestry informative markers in uterine neoplasms from women of European, African, and Asian ancestry.

Characterization of ancestry-linked peptide variants in disease-relevant patient tissues represents a foundational step to connect patient ancestry with disease pathogenesis. Nonsynonymous single-nucleotide polymorphisms encoding missense substitutions within tryptic peptides exhibiting high allele frequencies in European, African, and East Asian populations, termed peptide ancestry informative markers (pAIMs), were prioritized from 1000 genomes. In silico analysis identified that as few as 20 pAIMs can determine ancestry proportions similarly to >260K SNPs (R2 = 0.99). Multiplexed proteomic analysis of >100 human endometrial cancer cell lines and uterine leiomyoma tissues combined resulted in the quantitation of 62 pAIMs that correlate with patient race and genotype-confirmed ancestry. Candidates include a D451E substitution in GC vitamin D-binding protein previously associated with altered vitamin D levels in African and European populations. pAIMs will support generalized proteoancestry assessment as well as efforts investigating the impact of ancestry on the human proteome and how this relates to the pathogenesis of uterine neoplasms.

Personalized Single-Cell Proteogenomics to Distinguish Acute Myeloid Leukemia from Non-Malignant Clonal Hematopoiesis.

Genetic mutations associated with acute myeloid leukemia (AML) also occur in age-related clonal hematopoiesis, often in the same individual. This makes confident assignment of detected variants to malignancy challenging. The issue is particularly crucial for AML post-treatment measurable residual disease monitoring, where results can be discordant between genetic sequencing and flow cytometry. We show here, that it is possible to distinguish AML from clonal hematopoiesis and to resolve the immunophenotypic identity of clonal architecture. To achieve this, we first design patient-specific DNA probes based on patient's whole-genome sequencing, and then use them for patient-personalized single-cell DNA sequencing with simultaneous single-cell antibody-oligonucleotide sequencing. Examples illustrate AML arising from DNMT3A and TET2 mutated clones as well as independently. The ability to personalize single-cell proteogenomic assessment for individual patients based on leukemia-specific genomic features has implications for ongoing AML precision medicine efforts.

Tensin 1 (TNS1) is a modifier gene for low body mass index (BMI) in homozygous [F508del]CFTR patients.

Cystic fibrosis (CF) is a life-limiting autosomal recessive genetic disease caused by variants in the CFTR gene, most commonly by the [F508del] variant. Although CF is a classical Mendelian disease, genetic variants in several modifier genes have been associated with variation of the clinical phenotype for pulmonary and gastrointestinal function and urogenital development. We hypothesized that whole genome sequencing of a well-phenotyped CF populations might identify novel variants in known, or hitherto unknown, modifier genes. Whole genome sequencing was performed on the Illumina HiSeq X platform for 98 clinically diagnosed cystic fibrosis patient samples from the Adult CF Clinic at the University of California San Diego (UCSD). We compared protein-coding, non-silent variants genome wide between CFTR [F508del] homozygotes vs CFTR compound heterozygotes. Based on a single variant score test, we found 3 SNPs in common variants (MAF >5%) that occurred at significantly different rates between homozygous [F508del]CFTR and compound heterozygous [F508del]CFTR patients. The 3 SNPs were all located in one gene on chromosome 2: Tensin 1 (TNS1: rs3796028; rs2571445: and rs918949). We observed significantly lower BMIs in homozygous [F508del]CFTR patients who were also homozygous for Tensin 1 rs918949 (T/T) (p = 0.023) or rs2571445 (G/G) (p = 0.02) variants. The Tensin 1 gene is thus a potential modifier gene for low BMI in CF patients homozygous for the [F508del]CFTR variant.

Integrative Analysis of miRNAs Identifies Clinically Relevant Epithelial and Stromal Subtypes of Head and Neck Squamous Cell Carcinoma.

PURPOSE: The objective of this study is to characterize the role of miRNAs in the classification of head and neck squamous cell carcinoma (HNSCC). EXPERIMENTAL DESIGN: Here, we analyzed 562 HNSCC samples, 88 from a novel cohort and 474 from The Cancer Genome Atlas, using miRNA microarray and miRNA sequencing, respectively. Using an integrative correlations method followed by miRNA expression-based hierarchical clustering, we validated miRNA clusters across cohorts. Evaluation of clusters by logistic regression and gene ontology approaches revealed subtype-based clinical and biological characteristics. RESULTS: We identified two independently validated and statistically significant (P < 0.01) tumor subtypes and named them "epithelial" and "stromal" based on associations with functional target gene ontology relating to differing stages of epithelial cell differentiation. miRNA-based subtypes were correlated with individual gene expression targets based on miRNA seed sequences, as well as with miRNA families and clusters including the miR-17 and miR-200 families. These correlated genes defined pathways relevant to normal squamous cell function and pathophysiology. miRNA clusters statistically associated with differential mutation patterns including higher proportions of TP53 mutations in the stromal class and higher NSD1 and HRAS mutation frequencies in the epithelial class. miRNA classes correlated with previously reported gene expression subtypes, clinical characteristics, and clinical outcomes in a multivariate Cox proportional hazards model with stromal patients demonstrating worse prognoses (HR, 1.5646; P = 0.006). CONCLUSIONS: We report a reproducible classification of HNSCC based on miRNA that associates with known pathologically altered pathways and mutations of squamous tumors and is clinically relevant.

MutEnricher: a flexible toolset for somatic mutation enrichment analysis of tumor whole genomes.

BACKGROUND: Analysis of somatic mutations from tumor whole exomes has fueled discovery of novel cancer driver genes. However, ~ 98% of the genome is non-coding and includes regulatory elements whose normal cellular functions can be disrupted by mutation. Whole genome sequencing (WGS), on the other hand, allows for identification of non-coding somatic variation and expanded estimation of background mutation rates, yet fewer computational tools exist for specific interrogation of this space. RESULTS: We present MutEnricher, a flexible toolset for investigating somatic mutation enrichment in both coding and non-coding genomic regions from WGS data. MutEnricher contains two distinct modules for these purposes that provide customizable options for calculating sample- and feature-specific background mutation rates. Additionally, both MutEnricher modules calculate feature-level and local, or "hotspot," somatic mutation enrichment statistics. CONCLUSIONS: MutEnricher is a flexible software package for investigating somatic mutation enrichment that is implemented in Python, is freely available, can be efficiently parallelized, and is highly configurable to researcher's specific needs. MutEnricher is available online at https://github.com/asoltis/MutEnricher .

Molecular Analysis of Clinically Defined Subsets of High-Grade Serous Ovarian Cancer.

The diversity and heterogeneity within high-grade serous ovarian cancer (HGSC), which is the most lethal gynecologic malignancy, is not well understood. Here, we perform comprehensive multi-platform omics analyses, including integrated analysis, and immune monitoring on primary and metastatic sites from highly clinically annotated HGSC samples based on a laparoscopic triage algorithm from patients who underwent complete gross resection (R0) or received neoadjuvant chemotherapy (NACT) with excellent or poor response. We identify significant distinct molecular abnormalities and cellular changes and immune cell repertoire alterations between the groups, including a higher rate of NF1 copy number loss, and reduced chromothripsis-like patterns, higher levels of strong-binding neoantigens, and a higher number of infiltrated T cells in the R0 versus the NACT groups.

Neuroinflammatory astrocytes generated from cord blood-derived human induced pluripotent stem cells.

BACKGROUND: Astrocytes respond to central nervous system (CNS) injury and disease by transforming to a reactive astrogliosis cell state that can contribute to either CNS dysfunction or repair. Neuroinflammation is a powerful driver of a harmful A1 astrogliosis phenotype associated with in vitro neurotoxicity and histopathology in human neurodegenerative diseases. Here we report a protocol for the rapid development of a human cell culture model of neuroinflammatory astrogliosis using induced pluripotent stem cells (iPSCs). METHODS: Using RNA sequencing and in vitro cell assays, we measured transcriptional and cellular effects of chronic exposure of human iPSC-derived astrocytes to the cytokines TNFα (tumor necrosis factor alpha) or IL-1ß (interleukin-1 beta). RESULTS: We show TNFα and IL-1ß induce pro-inflammatory gene signatures but by widely different magnitudes. TNFα treatment results in 606 differential expressed genes, the suppression of glutamate-uptake, and increased phagocytic activity in astrocyte cultures. In contrast, IL-1ß effects are attenuated to 33 differential expressed genes and no significant effects on glutamate-uptake or increased phagocytic activity. CONCLUSION: Our approach demonstrates a rapid tool for modeling neuroinflammatory human astrocytic responses in nervous system trauma and disease. In particular, we reveal a model for robust TNFα-induced human astrogliosis suitable for the study of neurotoxic A1 astrocytes.